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Abstract 


Dimensional  analysis  and  dynamic  similarity  requirements  show  that  Reynolds,  Froude, 
and  Strouhal  scales  are  important  parameters  to  be  satisfied  for  free-running  submarine  model 
maneuvering  tests  in  order  that  they  correctly  represent  full  scale  behavior.  The  requirement  that 
all  three  scales  be  satisfied  simultaneously  presents  conflicting  demands  on  the  relationship 
between  the  velocity  U  and  the  characteristic  length  L  of  the  model.  The  inability  to  meet  these 
demands  presents  technical  challenges  to  properly  satisfy  dynamic  similarity  for  Radio- 
Controlled  Model  (RCM)  maneuvering  tests.  Instead,  we  take  a  first  principles  approach  by 
considering  the  forces  and  moments  acting  on  the  vehicle  and  seek  solutions  that  dynamically 
produce  the  same  scaled  forces  and  moments  on  the  model  as  exist  on  the  full  scale  vehicle.  In 
this  manner  we  dynamically  satisfy  the  above  three  scaling  parameters.  The  completion  of  the 
theory  will  lead  to  the  development  of  a  modified  test  procedure  for  and  subsequent  validation 
by  free-running  submarine  model  maneuvering  experiments. 

Summarizing,  we  will  develop  an  extended  dynamic  similarity  theory  that  can  be  applied 
to  the  unsteady  maneuvers  of  the  RCM.  We  will  apply  the  theory  practically  to  improve  the 
RCM  test  program  and  increase  correlation  with  full  scale  data.  This  will  provide  a  physical 
basis  for  deciding  which  improvements  are  better  made  in  simulation  and  also  enhance  the  model 
the  model  simulation  paradigm. 


Administrative  Information 

The  work  described  in  this  report  was  performed  by  the  Propulsion  and  Fluid  Systems 
Division  (Code  5400)  and  the  Maneuvering  and  Control  Division  (Code  5600)  of  the 
Hydromechanics  Department  at  the  Naval  Surface  Warfare  Center,  Carderock  Division 
(NSWCCD).  The  work  was  funded  by  the  Office  of  Naval  Research,  Code  331  as  part  of  the 
Scaling  Task  of  the  FY07  6.1  Turbulence  and  Stratified  Wakes  Program  (Program  Element 
0601 153N). 


Introduction 

The  Radio-Controlled  Model  (RCM)  is  a  semi-autonomous,  battery-powered,  free- 
running  1/1 6th  to  l/20th-scale  submarine  model  that  can  be  configured  as  a  model  of  any  of  the 
submarine  classes  in  the  fleet.  Because  it  is  free  running,  it  can  move  in  all  six  degrees  of 
freedom  and  provide  an  unsteady  maneuvering  capability.  The  RCM  can  perform  a  wide  range 
of  maneuvers  and  is  currently  the  best  predictor  of  full-scale  submarine  maneuvering 
performance.  It  is  an  integral  part  of  the  hydrodynamic  process  employed  by  the  Maneuvering 
and  Control  Division  to  develop  Submerged  Operating  Envelopes  and  other  hydrodynamic 
products  (Hess  et  al,  2004). 

When  RCM  maneuvering  simulation  data  are  compared  with  full-scale  data,  they  exhibit 
varying  degrees  of  Reynolds  scale  effect.  Forces  and  moments  acting  on  the  ship  hull, 
appendages  and  on  the  propulsor  are  known  to  experience  such  scale  effects.  With  a  two-order- 
of-magnitude  difference  in  Reynolds  number  between  model  and  full  scale,  one  should  not  be 
surprised  to  see  scale  effects  appearing  in  such  typical  RCM  maneuvers  as:  turns,  rudder  jam 
recoveries  and  stemplane  jam  recoveries. 
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To  see  how  the  Froude  scale  arises,  consider  a  deeply  submerged  submarine  undergoing  a 
maneuver  which  results  in  pitching  (and/or  rolling)  motion.  The  initial  phase  of  the  ship  motion 
response  is  characterized  by  the  excitation  moment  and  countered  by  the  restoring  moment.  The 
latter  is  a  product  of  ship  weight  and  a  moment  arm  produced  by  the  sine  component  of 
metacentric  height.  The  submarine  pitch  equation  of  motion,  when  made  dimensionless,  gives 
rise  to  a  dimensionless  quantity,  which  represents  the  ratio  of  the  hydrodynamic  exciting  moment 
to  the  hydrostatic  restoring  moment.  This  quantity  is  identical  to  the  Froude  number  as  it 
expresses  the  ratio  of  inertial  to  gravitational  effects  acting  on  the  submarine.  To  simulate  the 
restoring  moment  correctly  in  model  tests,  Froude  scaling  must  be  used,  and  this  dictates  the 
model  approach  speed  (Hess  et  al,  2004). 

As  a  free  running  model,  the  RCM  exhibits  unsteady  motion  during  maneuvering 
simulations.  Theodorsen  (1935)  shows  that  the  unsteady  dynamic  lift  and  phase  on  a  sinusoidal 
pitching  foil  oscillatory  motion  are  uniquely  characterized  by  the  reduced  frequency.  For  RCM 
maneuvers,  the  ship  motion  is  not  solely  oscillatory;  instead,  the  phase  and  unsteady  load  acting 
on  the  RCM  model  in  unsteady  maneuvering  are  characterized  by  the  Strouhal  number 
(Schlichting,  1979  and  Blake,  1984).  Strouhal  scaling  is  therefore  important  for  RCM 
maneuvering  simulation. 

In  addition  to  the  Strouhal  scale,  RCM  data  clearly  show  that  both  Froude  (  Fr )  and 
Reynolds  scales  ( Re )  are  needed  to  satisfy  the  dynamic  similarity  requirement,  yet  this  cannot 
be  accomplished  when  fresh  water  is  used  for  model  testing.  Fluid  material  parameters  such  as 
density,  p ,  dynamic  viscosity,  p ,  and  hence  kinematic  viscosity,  v ,  are  quite  similar  for  salt 
water  and  fresh  water.  Therefore,  when  matching  Reynolds  and  Froude  numbers  at  both  scales, 
one  must  chiefly  rely  on  the  variation  in  length  and  velocity  scales  because  the  fluid  material 
parameters  are  nearly  the  same.  Once  the  length  scale  is  set  by  the  chosen  model  scale  ratio, 
then  the  velocity  scale  is  used  to  effect  the  matching.  When  Froude  scaling  is  applied  to  set  the 
test  speed,  one  finds  that  the  Reynolds  number  is  80  to  100  times  greater  at  full  scale  than  at 
model  scale.  Consequently,  the  boundary  layer  thickness  and  boundary  layer  velocity  profile  are 
quite  different  between  full  scale  and  model  scale.  This  explains  the  unavoidable  appearance  of 
Reynolds  scale  effects  in  RCM  data. 

Extensive  test  programs  were  carried  out  in  the  1980’s  and  1990’s  to  improve  RCM 
simulation  test  results.  One  area  of  difficulty,  RCM  casualty  recovery  simulation,  showed 
insufficient  craft  pitch-up  motion  when  compared  with  full-scale  data.  A  method  to  enhance 
pitch  was  found,  and  it  involved  the  application  of  a  layer  of  heavy  grit  on  the  upper  surface  of 
the  hull.  Furthermore,  the  grit  appeared  to  have  no  untoward  effect  for  maneuvers  where  it  was 
not  normally  required.  This  grit  method  has  been  used  since  the  1990’s  for  all  RCM 
maneuvering  simulation  tests. 

Over  the  years,  however,  RCM  data  with  grit  has  displayed  mixed  results  with  good 
correlation  in  some  tests  and  poor  correlation  in  others.  The  mixed  results  with  grit  and  the 
fundamental  fact  that  Reynolds  numbers  cannot  be  matched  led  to  the  conclusion  that  the  test 
program  should  be  augmented  with  fixes  applied  in  simulation.  This  model  the  model  approach 
remains  the  current  standard  to  ensure  the  reliability  and  accuracy  of  the  existing  RCM  test 
method  to  simulate  full-scale  performance. 

While  the  RCM  remains  the  best  predictor  of  full-scale  submarine  maneuvering 
performance,  we  believe  that  improvements  to  the  test  program  can  be  made,  which  will  then 
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have  a  favorable  impact  for  modeling.  This  is  not  an  easy  task.  Clearly,  the  flow  field  around 
the  vehicle  is  quite  complex  during  maneuvering  and  particularly  during  casualty  recovery,  and  a 
reliable  numerical  model  requires  extensive  correlation  and  verification  with  model  and  full- 
scale  data.  We  observe  the  following  fundamental  issues  associated  with  the  existing  RCM  test 
method. 

•  The  physics  associated  with  the  use  of  a  grit  layer  to  improve  maneuvering  is  not 
understood. 

•  Froude  scaling  is  currently  used  to  set  the  model  approach  speed.  For  unsteady 
maneuvering,  the  requirement  should  be  that  the  Froude  number  is  matched  at  each 
instant  of  time  throughout  the  maneuver.  The  implication  is  that  the  ratio  of  the 
hydrodynamic  exciting  moment  to  the  hydrostatic  restoring  moment  may  not  remain 
correct  throughout  the  maneuver. 

•  Reynolds  scale  effects  cause  forces  and  moments  induced  by  full  scale  appendages  and 
the  hull  to  be  inadequately  simulated  in  RCM  tests. 

•  Strouhal  scale  which  is  an  important  parameter  influencing  unsteady  motions  is  not 
considered  in  RCM  tests. 

•  Existing  scaling  formulas  that  relate  model  data  to  full  scale  performance  are  not 
rigorously  derived. 

The  interrelated  consequences  of  each  of  these  issues  for  the  maneuvering  model  are  not 
understood.  Some  of  the  impacts  may  well  be  mitigated,  whereas  others  may  be  amplified.  This 
is  a  fundamental  question  that  needs  to  be  kept  in  mind.  As  the  theory  unfolds  in  the  Technical 
Approach  section,  the  challenge  will  be  to  apply  it  in  a  practical  manner  to  improve  not  only 
RCM  maneuvers,  but  also  the  resulting  model  the  model  simulation  process. 

Objectives 

Develop  an  extended  dynamic  similarity  theory  that  can  be  applied  to  the  unsteady 
maneuvers  of  the  RCM.  Apply  the  theory  practically  to  improve  the  RCM  test  program  and 
increase  correlation  with  full  scale  data.  Provide  a  physical  basis  for  improvements  that  are 
better  made  in  simulation  and  enhance  the  model  the  model  simulation  process. 

Technical  Approach 

Dimensional  analysis  and  dynamic  similarity  requirements  show  that  Reynolds,  Froude, 
and  Strouhal  scales  must  be  preserved  in  RCM  maneuvering  simulation  tests  (Baker,  1973; 
Gukhman,  1965).  Let  Re  denote  the  Reynolds  number,  which  is  defined  by  Re  =  UL/v ,  where 
U  and  L  denote  a  characteristic  speed  and  length,  respectively.  Define  Fr  as  the  Froude  number 
given  by  Fr  =  u/ -JgL  .  Finally,  St  denotes  the  Strouhal  number  represented  by 
St  =  {da/dt)L/U ,  where  da/dt  represents  a  characteristic  angular  velocity  (such  as  p,  q  or  r). 
The  Reynolds  number  describes  the  relative  importance  of  viscous  effects  to  inertial  effects,  the 
Froude  number  characterizes  the  ratio  of  hydrodynamic  to  hydrostatic  restoring  moments  acting 
on  the  model,  and  the  Strouhal  number  is  a  fundamental  measure  of  the  unsteadiness  of  the 
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motion.  The  subscripts  s  and  m  will  denote  the  full  scale  ship  and  the  model,  respectively.  We 
then  wish  to  obtain: 


(UL/  v)m  =  (UL/V)s  Re  scale 

Fr  scale  .  (1) 

[(da/dt)L/U]m  =  [(da/dt)L/u]s  St  scale 

These  equations  show  that  since  U  and  L  appear  both  in  the  numerators  and  in  the  denominators 
of  these  quantities,  they  cannot  be  matched  simultaneously.  We  have  developed  a  new  approach 
to  address  this  technical  challenge. 


Basic  Concepts 

Ship  motion  is  a  response  to  forces  and  moments  acting  on  a  craft.  A  typical  submarine 
configuration  showing  rudders,  stemplanes,  sail  and  the  general  shape  of  the  hull  is  depicted  in 
Fig.  1.  The  classic  approach  of  satisfying  geometric  and  kinematic  similarities  in  model  tests  is 
intended  to  produce  similar  forces  and  moments  at  model  and  full  scale  to  generate  the  same  ship 
motions.  Generating  proper  forces  and  moments  acting  on  the  model  is  the  first  priority  in 
conducting  model  tests. 


Figure  1.  A  Typical  Submarine  Configuration 


The  forces,  moments  and  the  corresponding  ship  motion  can  be  expressed  in  terms  of 
three  linear  and  three  angular  differential  equations  similar  in  form  to  that  which  describes  a 
spring-mass  system.  Furthermore,  Lloyd  (1989)  has  shown  that  for  conventional  surface  ships  in 
waves  the  product  moments  of  inertia  are  usually  small  and  negligible.  Let  Fi  and  Mt  denote 

force  and  moment  components,  which  will  be  normalized  by  \/2pU2S  and  l/ 2 pU2S  L , 
respectively,  where  we  use  the  symbols  L  and  S  to  denote  ship  length  and  hull  surface  area.  Let 
m0  and  70  denote  ship  mass  and  moment  of  inertia,  and  m,  and  7,  denote  virtual  mass  and 

moment  of  inertia.  The  quantities  d2xjdt 2  and  d2ajdt2  will  denote  linear  and  angular 
acceleration  components.  The  index  variables  /  =  1,2,3  and  i  =  4,5,6  will  represent  surge,  sway 

and  heave  for  linear  motions  and  roll,  pitch  and  yaw  for  angular  motions,  respectively.  We  then 
obtain: 
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i  =  1,2,3  and 
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i  =  4,5,6 


If  normalized  forces  and  moments  on  the  left-hand-sides  (LHS)  of  Eqs.  2  are  the  same  for 
both  the  model  and  the  ship,  namely  LHSm  =  LHSs ,  then 
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Equality  of  the  forces  and  moments  will  then  result  in  a  vehicle  response,  which  must  be 
the  same  at  model  and  full  scale;  namely,  the  inertia  terms  on  the  right-hand-sides  (RHS)  of 
Eqs.  2  will  be  the  same  and  be  given  by 
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i  =  1,2,3  and 
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From  Eqs.  4,  scaling  formulas  or  transfer  functions  can  be  derived  to  relate  model  data  to 
full  scale  ship  motion.  The  challenge  will  be  to  develop  an  experimental  procedure  for 
submarine  model  testing  that  can  produce  the  same  normalized  forces  and  moments  acting  at 
both  scales  as  specified  in  Eqs.  3. 


Mathematical  Formulation 

First  consider  the  angular  motions.  The  linear  motions  will  be  treated  later.  Let  a 
denote  an  Euler  angle  component,  where  a,  =  (<p,0,y/)  and  (p ,  0  and  y/  represent  roll,  pitch, 

and  yaw  angles,  respectively.  The  ship  motion  can  be  generalized  in  the  form  of  a  differential 
equation  for  a  spring-mass  system  for  each  angle  component  as  follows  (refer  to  Fig.  1). 


Me+M](a)  +  M2 


(da > 

=  M, 

( d2a ) 

<  dt  ; 

J 

[dt2] 

=  (/„+/,)■ 


d2a 
dt 2 


(5) 
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The  term  Me  refers  to  excitation  moments  from  movable  appendages  such  as  rudders, 
stemplanes  and  a  propulsor.  The  terms  M2(da/dt)  and  Mi(d2a/dt2)  denote  the 

moments  that  are  functions  of  flow  angle,  angular  velocity,  and  angular  acceleration, 
respectively.  Equation  5  applies  both  at  full  scale  and  model  scale,  namely 


M,+M,(a)+M,W 

J  5 

(t  r\^2« 

M,  +  */,«/)+ M,^] 

-  m 

/  \d~ct 

[(°  ')<*’] 

(6) 


The  first  task  in  this  study  is  to  identify  a  test  method  such  that  forces  and  moments 
acting  on  the  full  scale  vehicle  are  correctly  represented  at  model  scale,  namely  to  achieve 


(  da} 

Me+M,(a)  +  M2  — 
\dt  ) 

A/e+M,(«)  +  M2^ 

pU2SL 

m 

—  pU2SL 

2 

The  next  task  will  be  to  derive  the  scaling  formulas  to  relate  measured  RCM  data  to  full 
scale  ship  motions.  This  task  is  performed  by  equating  the  vehicle  response  on  the  right-hand- 
sides  of  Eqs.  6 
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(8) 


Test  Procedure  to  Satisfy  Dynamic  Similarity  Requirements 

We  will  assume  that  forces  and  moments  acting  on  a  full  scale  vehicle  can  be  simulated 
with  RCM  tests  if  Reynolds,  Froude,  and  Strouhal  numbers  are  preserved  at  each  scale. 
Simulation  of  a  ship  in  a  planes-fixed  turn  will  be  used  as  an  example  to  demonstrate  the 
concepts  and  to  develop  a  new  RCM  test  procedure.  A  similar  approach  can  be  applied  to  rudder 
or  stemplane  jam  recovery  simulations.  Note  that  a  planes-fixed  turn  is  one  for  which  no 
automatic  controller  is  used  with  the  stemplanes  to  maintain  a  constant  depth  during  the  turn; 
instead,  the  stemplanes  are  maintained  in  a  fixed  position,  and  the  depth  is  allowed  to  vary. 

A  sketch  of  craft  trajectory  and  the  corresponding  response  of  the  vehicle  as  described  by 
time  histories  of  ship  speed,  angular  velocity,  and  accelerations  in  response  to  a  rudder  angle 
deflection  are  shown  in  Fig.  2  (Lewis,  1989).  The  terms  v,  /?  and  r  denote  lateral  velocity,  drift 
angle  and  yaw  angular  velocity  (turning  rate).  The  terms  v  and  r  represent  accelerations.  The 
rudder  is  deflected  with  a  constant  deflection  rate  of  ddr  /dt  to  a  specified  angle  Sr  as  shown  in 
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Figure  2.  Turning  Maneuver  Description  and  Phases 
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Fig.  2.  The  craft  motion  is  divided  into  three  phases  with  the  start  of  the  rudder  deflection 
representing  the  time  origin,  /  =  0 . 

The  first  phase  is  dominated  by  the  rudder  deflection,  and  the  ship  response  is  primarily 
inertial.  The  second  phase  is  associated  with  ship  turning  represented  by  the  appearance  of  side 
velocity  v  and  turning  rate  r.  The  quasi-steady  turn,  as  demonstrated  by  the  disappearance  of 
inertial  terms,  signifies  the  third  phase.  To  show  a  turning  example,  the  symbols  used  in  Eq.  6 
are  revised  as 


Me  +  M,(v)  +  M2(r) 

Me  +  M,(v)  +  M2(r) 

~pU2SL 

m 

pU2SL 

(7a) 


where  v  is  the  lateral  velocity  at  the  center  of  gravity  (CG)  and  is  obtained  from  v  =  —U  sin  /? . 

Motivation  for  Reynolds  Scale  Simulation  in  RCM  Tests 

The  requirement  to  match  Reynolds  numbers  results  from  a  consideration  of 
Conservation  of  Momentum  for  the  fluid  in  which  the  vehicle  is  moving.  This  physical  law  is 
expressed  in  the  Navier-Stokes  equations,  and  when  these  are  made  dimensionless,  the  Reynolds 
number  arises  as  the  key  dimensionless  group  expressing  the  ratio  of  inertial  to  viscous  forces. 
By  matching  this  dimensionless  group  at  both  model  and  full  scale,  the  governing  equations  for 
the  fluid  will  then  be  identical  and  will  result  in  similar  solutions.  As  discussed  above  in  Eq.  1, 
we  must  also  match  Froude  and  Strouhal  numbers.  As  we  proceed  we  will  find  that  it  is  possible 
to  match  the  latter  two  groups,  but  we  cannot  also  match  the  Reynolds  number  simultaneously. 
Therefore,  the  two-order-of-magnitude  difference  in  Reynolds  number  between  model  and  full 
scale  implies  Reynolds  scaling  effects  will  appear  at  model  scale  for  those  motions  where 
viscous  effects  are  important.  The  impact  of  these  effects  will  be  to  alter  the  force  and  moment 
distribution  acting  on  the  vehicle  to  various  degrees.  Although  we  cannot  match  Reynolds 
number,  if  we  can  assess  how  these  force  and  moment  changes  are  created  by  the  appendages 
and  the  hull,  then  we  can  correct  them.  We  will  refer  to  this  correction  approach  as  dynamically 
satisfying  Reynolds  number.  The  next  few  subsections  will  evaluate  the  needed  corrections  on 
various  appendages  and  the  hull  for  a  submarine  turning  maneuver. 

Reynolds  Scale  Correction  on  Rudder  Performance  -  First  Phase 

As  noted  in  Fig.  2,  the  first  phase  of  a  turning  maneuver  is  that  short  period  of  time  when 
the  vehicle  motion  is  still  governed  by  inertia  and  corresponds  to  v  »  0 ,  r  ~  0,  v  *  0  and  r  *  0 . 
We  can  simplify  Eq.  7a,  to  obtain 


Me 

K 

—  pU2SL 

12  J 

m 

X-pU2SL 

(9) 
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Equation  9  requires  that  the  dimensionless  moments  produced  by  rudders  to  turn  the  vehicles 
about  their  CG  must  be  equal  if  the  first  phase  of  the  turning  maneuver  is  to  be  properly 
simulated  in  RCM  tests. 

As  can  be  seen  in  Fig.  1,  rudders  are  located  at  the  aft  end  of  the  vehicle  and  are 
immersed  in  the  boundary  layer  developed  along  the  ship’s  hull.  The  ratio  of  the  boundary  layer 
thickness  to  the  ship  length  is  much  greater  for  the  model  than  for  the  full  scale  vehicle  as  shown 
by  the  typical  profiles  given  in  Fig.  3.  The  model  rudder  placed  in  the  thicker  boundary  layer 
will  develop  a  smaller  force  and  moment  than  a  full  scale  rudder.  This  Reynolds  scale  effect  is 
associated  with  the  boundary  layer  velocity  deficit. 


Figure  3.  Boundary  Layer  Velocity  Profiles  at  the  Rudder  Plane 

Another  scale  effect  acting  on  the  rudder  is  related  to  the  boundary  layer  displacement 
thickness.  The  Reynolds  numbers  based  on  mean  rudder  chord  are  around  Re  =  0.5  x  106  for  the 
model  and  Re  =  40  x  106  for  a  full  scale  rudder.  Measurements  on  several  NACA  airfoils,  shown 
in  Fig.  4,  exhibit  significant  scale  effects  on  lift  coefficients  especially  at  low  Reynolds  numbers 
of  less  than  1  x  106  (Loftin  and  Smith,  1949;  Loftin  and  Bursnall,  1950;  and  Whicker  and 
Fehlner,  1958).  The  lift  coefficients  are  seen  to  increase  with  Reynolds  number.  This  trend  is 
also  observed  theoretically  by  Spence  (1954). 

Consider  a  simple  rectangular  foil  of  chord  length  c  and  span  s  located  in  a  non-uniform 
flow  u(y)  as  shown  in  Fig.  5.  The  section  lift  A£  on  a  segment  Ay  due  to  an  angle  of  attack  8r 
is  given  by 


A  £  =  Q 


cAy  =  d£±Sr 

dSr  r 


cAy . 


Taking  the  limit  as  Aj>  — >  0 ,  we  obtain 


_  f  SCL 


dS, 


^/>«2O0 


cdy  =  ^ ocSr  J ^r(y)  u\y)dy , 


(10) 


(11) 


9 


Section  lUrt-ourve  slope  per  degree,  doj/da, 


Figure  4.  Variation  of  Lift  Curve  Slope  of  17  Airfoils 
of  Smooth  Surfaces  with  Reynolds  Numbers 
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SttotiloQ  lift-curve  slope  par  degree,  doj/dOg 
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Figure  4  (Cont.).  Variation  of  Lift  Curve  Slope  of  17  Airfoils 
of  Roughened  Surfaces  with  Reynolds  Numbers 
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where  dCL/dSr  (y)  denotes  the  lift  curve  slope,  which  varies  along  the  span  as  a  result  of  the 
boundary  layer  velocity  gradient.  However,  the  weak  variation  in  dCL/dSr  (y)  in  the  boundary 
layer  at  model  scale  is  much  smaller  than  the  variation  resulting  from  a  two-order-of-magnitude 
Reynolds  number  difference  between  model  and  full  scale.  Hence,  we  will  replace  the  weak 

boundary  layer  variation  with  a  mean  value,  dCL/dSr ,  and  remove  it  from  the  integral. 


Figure  5.  Boundary  Layer  Effect  on  Lifting  Surfaces 


We  obtain 


t=\pcS'  <12) 

Shift  now  to  dimensionless  variables  u  =  u/U  and  y'  =  y/s  to  get 

l  =  ^PcSrTtUls  fn  (13) 

2  oor  J0 

Let  the  rudder  area  be  denoted  by  Ar  =  cs  and  divide  both  sides  by  hull  surface  area  S  to  get  a 
normalized  lift  given  by 

£ 

\PU‘S 

Use  Xr  to  denote  the  distance  between  the  rudder  and  the  center  of  gravity  of  the  vessel. 
Multiply  both  sides  by  Xr/L  to  get  a  normalized  turning  moment  in  the  form 


— f  u'2(y')dy' . 
e  r  nx  Jo  v' 


dS. 
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(15) 


where  the  subscripts  are  used  to  represent  the  equations  evaluated  at  full  and  model  scales. 

By  now  substituting  the  above  equations  into  Eq.  9  and  noting  that  [XjL]  t  =  \XjL]m  we 
obtain 


(16) 


where  £  denotes  the  Reynolds  scale  effect  on  rudder  lift  as  a  result  of  boundary  layer 
displacement  thickness,  and  7  the  scale  effect  due  to  the  difference  in  the  hull  boundary  layers  at 
each  scale.  Summarizing,  we  have 


W.  ■[*,]. 
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dS, 


dCL 
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£  ",2(y'W 


£  u'\y'W 


dCL 
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SCL 
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(17) 


A  Reynolds-averaged  Navier-Stokes  (RANS)  code  was  used  to  calculate  boundary  layer 
velocity  profiles  on  three  vessels.  The  numerical  results  for  7  are  shown  in  Table  1  and  are 
found  to  be  sensitive  to  hull  shape  with  values  ranging  from  1.06  to  1.16. 
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Table  1.  Values  for  ?7  =  ^{0  w'2(/)^'j  At  u'\y')dy' 


Top  Rudder 

Bottom  Rudder 

Ship  I 

1.06 

1.09 

Ship  II 

1.12 

1.16 

Ship  III 

1.07 

1.08 

The  effect  of  Reynolds  scale  on  lift  curve  slope  £  has  been  calculated  theoretically  by  Spence  on 
a  two  dimensional  Joukowski  airfoil  at  a  6°  attack  angle  and  is  shown  in  Fig.  6.  The  value  of  £ 
is  estimated  to  be  around  £  *  1.13  for  a  model  Reynolds  number  of  Re  =  0.5  x  106  and  a  full 
scale  value  of  Re  =  40  x  106 .  Thus,  a  significant  viscous  effect  on  rudder  force  is  noted.  Two 
methods  to  satisfy  Eq.  16  are  presented. 


Figure  6.  Reynolds  Number  Effect  on  Lift  Slope 


Method  1 

Modify  the  rudder  deflection  angle  in  the  RCM  test.  For  this  case,  the  model  rudder  is 
constructed  to  be  geometrically  similar  to  that  for  full  scale,  namely  [^r/^]5  —  [v4r/Sjw ,  and 

Eq.  16  becomes 

foL  =  [A-L  £  »7»  where  4  >  1  a°d  t]  >  1 .  (18) 
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As  a  numerical  example,  consider  the  case  of  a  full  scale  vessel  conducting  a  turning  maneuver 
with  a  rudder  deflection  of  [Sr ] s  =  25° .  To  compensate  for  the  Reynolds  scale  effect  on  rudder 

force  and  turning  moment,  the  model  rudder  must  be  deflected  to 

l^r]m  =[^r]s  £7  =  (25Xl-13Xl-06)=30°  for  one  hull  shape,  and 
[^r  ]  m  =  [^r  L  £  7  =  (25Xl  •  1 3Xl  •  1 6)  =  33°  for  another  hull  shape 

Note  carefully  that  the  values  of  £  and  77  used  in  this  example  are  just  for  demonstration.  More 

accurate  values  will  need  to  be  obtained  from  a  correlation  study  of  foil  scale  and  RCM 
maneuvering  data. 


Method  2 

Modify  the  rudder  size  for  the  RCM  model.  In  this  approach,  the  rudder  deflection 
angles  for  execution  of  the  turn  will  be  the  same  for  RCM  tests  and  for  foil  scale  trials.  That  is, 
we  will  set  [Sr]m  =  [<5'r]j ,  and  Eq.  16  becomes 


A 


£t]  =  l  .or 


kl 

rar  ic 


#7= 


one  hull  shape 
another  hull  shape 


(20) 


where  A  =  LS/Lm  is  the  length  scale  ratio  between  foil  and  model  scales.  This  implies  that  the 

scale  effect  on  rudder  force  and  moment  can  be  mitigated  if  the  rudder  on  the  RCM  model  is 
made  20-3 1  %  greater  in  area  than  a  geometrically  similar  rudder. 


Second  Phase  of  a  Turning  Maneuver 

The  second  phase  of  a  turning  maneuver  corresponds  to  conditions  where  the  vehicle 
begins  to  respond  to  the  rudder  deflection  with  lateral  and  yawing  motions  such  that 
v  *  0,  r  *  0,  dv/dt  *  0  and  dr/dt  *  o .  Under  these  conditions,  we  return  to  Eq.  7a  to  establish 
the  dynamic  scaling  requirement: 


Me  +  M,(v)  +  M2(r) 

Me  +  A/,(v)  +  M2(r) 

- pU2SL 

2  j 

m 

^ pU2SL 

(7a) 


The  presence  of  nonzero  v  and  r  terms  will  produce  varying  flow  attack  angles  along  the 
ship  hull,  sail  and  control  appendages,  such  as  rudders  and  stemplanes. 
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Reynolds  Scale  Correction  on  Sail  Performance 

At  model  scale,  the  sail  typically  operates  at  a  Reynolds  number  of  Rem  =  1 .5  x  106 .  Let 
u,  v,  and  w  denote  velocity  components  at  the  ship  center  of  gravity  (CG).  The  angular  velocities 
for  roll,  pitch,  and  yaw  motions  will  be  given  by  p,  q,  and  r.  Let  Xsai,  be  the  distance  between 

the  CG  of  the  ship  and  the  sail,  such  that  it  represents  the  moment  arm  for  turning.  ZsaU ,  a 
positive-definite  scalar,  will  denote  the  vertical  distance  between  the  location  of  the  resultant 
force  produced  by  the  sail  and  the  x-axis  of  the  ship,  thereby  denoting  a  moment  arm  for  roll. 

The  side  slip  angle  /?  at  the  sail  can  be  calculated  by 


vsail=v  +  Pzsail+rXsail 


(21) 


Dynamic  lift  will  be  produced  by  the  sail  due  to  the  angle  of  attack,  P .  A  sail  is 
typically  located  forward  of  the  CG  of  the  ship  where  the  boundary  layer  that  develops  along  the 
ship’s  hull  is  still  relatively  thin.  Using  a  typical  l/7th  power  law  turbulent  boundary  layer 
velocity  profile,  a  value  of  7  » 1 .02  may  be  estimated  on  the  ONR  Body  1  sail  using  available 
geometrical  information,  (Faller  and  Merrill,  2001).  A  value  of  £  ~  1.07  is  estimated  by 
applying  Spence’s  data  provided  in  Fig.  6  for  Rem  =  1 .5  x  106 .  Using  these  two  values  for  £  and 
7  and  the  method  described  in  Eq.  20,  we  can  correct  the  Reynolds  scale  effect  on  the  sail 


kL  _  [s]. 
kL =  W. 


£7  =  1 .09/ A2 . 


(22) 


Equation  22  indicates  that  the  scale  effect  on  the  force  and  moment  produced  by  the  sail 
on  the  RCM  model  can  be  mitigated  by  making  the  model  sail  9  %  greater  than  a  geometrically 
similar  sail. 


Reynolds  Scale  Correction  on  Ship  Hull  Performance 

At  model  scale,  the  ship’s  hull  typically  operates  at  a  Reynolds  number  of 
Re„,  =  1 .5  x  107 .  Steady  flows  over  axisymmetric  bodies  at  an  angle  of  attack  have  received 

considerable  attention  both  theoretically  and  experimentally  in  the  areas  of  three  dimensional 
flow  separation,  the  formation  of  leeward  vortices  and  the  unsteadiness  around  separation  zones 
(Bridges,  2003;  Fu  et  al,  1994;  Han  and  Patel,  1979).  A  6:1  prolate  spheroid  undergoing  time 
dependent  maneuvers  was  tested  by  Hoang,  et  al.  (1994)  using  their  Dynamic  Plunge-Pitch-Roll 
(DyPPiR)  facility  at  a  Reynolds  number  of  Rem  =  4.2  x  106 .  Their  data  showed  that  unsteady 

flow  separation  phenomena  were  significantly  different  from  those  associated  with  steady 
separation.  Furthermore,  airfoils  at  an  angle  of  attack  in  unsteady  flow  exhibit  significantly 
delayed  flow  separation  (Shen  and  Fuhs,  1999).  Lloyd  (1989)  observed  on  surface  ships  that 
viscous  forces  play  a  less  important  role  in  ship  motion  dynamics.  Two  turning  cases  are 
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considered.  One  corresponds  to  a  gentle  turn  with  small  v  and  r  terms,  and  the  other  corresponds 
to  a  tight  turn  with  large  v  and  r  terms. 

Gentle  Turn 

There  is  no  scale  effect  on  the  ship’s  hull  due  to  the  boundary  layer  velocity  deficit  term, 
namely  7  =  1.  It  is  further  noted  from  unsteady  experimental  data  that  the  flow  angles  associated 
with  flow  separation  observed  in  steady  tests  are  greatly  delayed.  For  a  gentle  turn,  drift  angle  at 
all  locations  along  the  hull  is  small,  and  cross  flow  may  not  play  a  significant  role  in  producing 
force  and  moment.  Based  on  the  theories  of  Spence  (1954),  Howarth  (1935)  and  Bridges,  et  al. 
(2003),  £  «  1.01  to  1.02  is  estimated  for  the  ship  hull  at  Rem  =  1.5  x  107 .  Using  these  two  values 
for  £  and  7  and  the  method  described  in  Eq.  20,  the  Reynolds  scale  effect  on  the  hull  may  be 
estimated  to  be 


\A  1 

A2T^  =  ^7  =  10lto,.°2.  (23) 

\Ahuu\s 

The  scale  effect  on  ship  hull  is  small.  This  is  in  agreement  with  general  observations  in 
straight  line  and  rotating  arm  data.  If  model  experiments  are  performed  at  Reynolds  numbers 
greater  than  10  to  15  million,  then  scale  effects  between  model  and  full  scale  appear  to  be 
negligible  for  the  purpose  of  making  stability  and  control  predictions. 

Tight  Turn 

During  a  tight  turn,  accompanied  by  larger  v  and  r  terms,  large  flow  angles  as  much  as  20 
to  30  degrees  may  be  encountered  at  the  ship  stem.  Therefore,  cross  flow  separation  may  be  an 
issue.  The  degree  to  which  scale  effects  play  a  role  in  the  cross  flow  separation  for  RCM 
maneuvering  tests  at  Rem  =  1.5  x  107  is  not  clear.  This  is  an  area  which  warrants  further  study. 

Surface  roughness  (grit)  has  been  applied  to  the  RCM  hull  and  has  demonstrated  considerable 
influence  on  turning  maneuvers.  If  scale  effects  on  cross  flow  separation  should  prove  to  be 
large  in  RCM  experiments,  then  a  closer  examination  of  distributed  hull  roughness  should  be 
carried  out.  Another  possibility  is  the  application  of  a  thin  wire  along  the  side  of  the  ship  hull  to 
trigger  a  supercritical  flow  pattern  as  demonstrated  in  a  sphere  experiment  (Schlichting,  1979). 
The  resulting  transition  to  turbulent  cross  flow  can  greatly  modify  cross  flow  drag. 

Reynolds  Scale  Correction  on  Rudder  Performance  -  Second  Phase 

Equation  7a  describes  the  moments  to  be  matched  between  model  and  full  scales.  The 
simpler  Eq.  9  was  obtained  by  removing  terms  that  are  negligible.  The  remaining  term,  the 
exciting  moment,  results  from  the  rudder  deflection.  During  the  second  phase  of  the  turn, 
nonzero  lateral  and  yawing  motions  modify  the  moment  produced  by  the  rudder  deflection  such 
that  it  may  be  said  to  originate  from  a  smaller  effective  rudder  angle.  Further,  this  effective 
rudder  angle  changes  continuously  in  time  (grows  smaller)  while  the  ship  is  in  a  turn  as  v  and  r 
grow  in  magnitude.  Let  Sr  denote  the  geometric  rudder  angle  with  respect  to  the  ship  centerline. 
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The  effective  rudder  angle  Sreff  is  defined  by  the  difference  of  the  geometric  rudder  angle  and  a 
correction  term  which  represents  the  instantaneous  flow  angle  acting  on  the  rudder 


Sreff  =  Sr-c(v,r) 


=  5  -  tan 


-i 


-\-rX, 


(24) 


where  Xr ,  a  positive-definite  scalar,  is  the  distance  from  the  CG  aft  to  the  rudder.  Following 
Eq.  18  for  geometrically  similar  rudders,  we  now  use  the  effective  rudder  angle  to  get 


[Sreff\m=[Sreff\jn 


[8r-c\m=[8r-c\'tl 

[sr L  -  [4  =  teL  -  [4  k  ?  but  W-  =  [4 =  c  • 

[Sr]m=[Sr]jrj-c{4rj-\) 


foL  =fc],£  7 -tan"1 


-v-rX. 


(Zn-i) 


(25) 


The  last  of  Eqs.  25  gives  the  angle  to  which  the  rudder  should  be  set  at  model  scale  to  produce 
the  same  exciting  moment  as  that  obtained  at  full  scale.  The  correction  for  scale  effects  given  by 
Eq.  18  for  phase  1  and  by  Eq.  25  for  phase  2  should  be  applied  at  each  instant  of  time.  However, 
during  phase  2,  as  v  and  r  grow  in  magnitude,  the  effective  rudder  angle  decreases  as  does  the 
exciting  moment  implying  that  the  need  for  scale  effect  correction  diminishes.  Thus,  we  propose 
a  rudder  angle  correction  as  shown  in  Fig.  7,  which  is  larger  during  Phase  1  and  the  early  part  of 
Phase  2,  but  decreases  in  magnitude  as  the  turn  proceeds.  The  details  of  the  shape  of  the 
correction  can  be  refined  after  conducting  model  tests. 


Figure  7.  Proposed  Rudder  Deflection  in  Model 
Tests  for  Reynolds  Scale  Correction 
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Table  2  summarizes  the  estimated  magnitudes  of  Reynolds  scale  effect  corrections  on  the 
exciting  moments  produced  by  rudders,  stemplanes,  the  sail,  and  the  ship  hull. 

Table  2.  Reynolds  Scale  Effects 


Item 

$ 

7 

Sn 

Rudder 

1.13 

1.06-1.16 

1.20-1.31 

Stemplane 

1.13 

1.06-1.16 

1.20-1.31 

Sail 

1.07 

1.02 

1.09 

Hull 

1.01-1.02 

1.00 

1.01-1.02 

Motivation  for  Froude  Scale  Simulation  in  RCM  Tests 

For  submarines,  the  size  of  the  top  rudder  is  typically  greater  than  the  bottom  rudder. 
Both  rudders  are  linked  such  that  they  respond  with  the  same  rudder  angle  when  executing  a 
turn.  When  a  mm  is  initiated,  a  roll  moment  is  produced  from  the  different  size  rudders,  and  the 
ship  will  experience  roll  motion  in  addition  to  yawing.  A  hydrostatic  restoring  moment  is  then 
induced  resulting  from  the  couple  produced  by  weight  acting  at  the  CG  and  buoyancy  acting  at 
the  center  of  buoyancy  (CB).  Buoyancy  is  expressed  by  B  =  pVg ,  where  V  is  the  volume  of  the 
ship.  For  a  submarine,  the  CB  is  located  a  distance  ZBG  directly  above  the  CG,  and  this 
metacentric  height  (reduced  by  the  sine  of  the  roll  angle)  serves  as  the  moment  arm  for  the 
restoring  moment.  Let  cp  denote  the  ship  roll  angle.  The  similarity  requirement  to  ensure  the 
same  dimensionless  restoring  moments  during  the  mm  can  be  expressed  by 


pVgZBG  sin  <p 

pVg  ZBG  sin  (p 

\/2pU2SL 

m 

\/2pU2SL 

(26) 


Eq.  26  will  be  satisfied  if 


S^BG 

S^BG 

U2  \ 

m 

L  u2  J 

and 


M«  =  i<P\s  ■ 


(27) 


The  first  condition  of  Eq.  27  states  that  the  model  must  be  tested  at  Froude  scale.  A  test 
method  to  ensure  that  Froude  scale  is  preserved  during  RCM  testing  will  be  discussed  later.  The 
second  condition  of  Eq.  27  requires  that  the  angle  experienced  during  roll  motion  at  model  scale 
must  be  the  same  as  that  measured  at  full  scale.  The  conditions  needed  to  satisfy  these 
requirements  during  RCM  testing  will  be  examined. 

Motivation  for  Strouhal  Scale  Simulation  in  RCM  Tests 

A  ship  in  a  mm  develops  angular  velocities  in  roll,  pitch  and  yaw  associated  with  the 
M2(r)  term.  Unsteady  lift  and  moment  are  produced  as  a  result.  The  dynamic  lift  on  an  airfoil 
in  a  sinusoidal  pitching  motion  has  been  theoretically  investigated  in  an  inviscid  fluid  by 
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Theodorsen  (1935).  His  theory  shows  that  the  dynamic  lift  and  the  phase  angle  with  respect  to 
the  instantaneous  foil  angle  is  uniquely  defined  by  the  reduced  frequency  K  defined  by 


(28) 


where  to ,  C,  and  U  denote  the  oscillation  frequency  of  the  attack  angle  variation,  chord  length, 
and  free  stream  speed,  respectively. 

During  RCM  maneuvering  tests,  the  motion  in  the  second  phase  of  the  turning  maneuver 
is  unsteady.  The  angular  velocity  da/dt  is  not  constant,  but  is  changing  continuously  with  time. 
Experience  shows  that  the  unsteady  forces  and  moments  can  be  characterized  by  a  Strouhal 
number  St  as  follows  (Schlichting,  1979  and  Blake,  1984), 

St  =  —~,  (29) 

dt  U 


Matching  the  Strouhal  scale  requires  that 


da  L 

da  L 

~dtU_ 

m 

~dtU. 

(30) 


Equation  30  is  an  important  equation  to  be  satisfied  in  model  tests,  if  the  amplitude  and  phase  of 
the  unsteady  forces  and  moments  are  to  be  correctly  simulated  at  model  scale.  To  satisfy  Eq.  30, 
we  will  make  a  prudent  choice  that 


da 

da 

_  dt . 

m 

,dt  _ 

Mm. 
•  U. 


(31) 


Substituting  the  right-hand  side  of  Eq.  3 1  into  the  left-hand  side  of  Eq.  30  will  then  dictate  that: 


(32) 


Therefore,  the  angular  velocity  to  be  experienced  by  the  RCM  model  must  be  maintained  at  an 
angular  velocity  experienced  by  the  ship  multiplied  by  a  factor  of  UjUm  .  We  will  see  in  a  later 
section  that  Eq.  32  can  be  satisfied  by  matching  the  Froude  number  and  that  Eq.  3 1  will  then 
follow  as  a  result  of  this  matching. 

These  are  the  requirements  that  ensure  that  unsteady  motions  are  correctly  simulated  in 
RCM  tests.  Methods  to  satisfy  these  two  requirements  will  be  discussed  in  a  later  section. 

Crashback  in  Sternplane  (Rudder)  Jam  Recovery  Simulation 

During  the  process  of  recovery  from  a  sternplane  jam,  a  large  rudder  deflection  is 
typically  executed  to  turn  the  ship  sideways  to  produce  a  large  drift  angle  and  large  form  drag  to 
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slow  the  ship.  This  is  often  followed  by  a  rapid  reversal  of  direction  of  propeller  rotation  (RPM 
reversal)  to  produce  a  large  negative  thrust  to  decelerate  and  stop  the  ship.  This  type  of 
maneuver  where  the  propeller  is  in  reverse  but  the  vehicle  still  has  forward  velocity  is  termed 
crashback.  Thus,  a  major  distinction  between  a  conventional  turning  maneuver  and  a  crashback 
maneuver  is  the  propeller  RPM  reversal.  As  an  example  we  will  use  a  model  test,  which  is 
conducted  for  the  purpose  of  correlating  with  known  full  scale  crashback  maneuvers  conducted 
during  trials  of  the  full  scale  vessel.  Therefore,  we  assume  that  time  histories  of  full  scale  ship 
speed,  RPM,  velocity  components,  attitude  and  trajectory  are  known.  Generally,  after  extensive 
earlier  predictions  prior  to  full  scale  trials  are  provided  to  the  fleet,  a  series  of  model  scale 
casualty  recovery  maneuvers  are  carried  out  after  full  scale  trials  to  validate  the  reliability  of  the 
RCM  simulation  capability.  The  procedure  discussed  here  makes  use  of  the  latter  as  an  example, 
but  the  method  can  be  reversed  so  that  RCM  tests  can  be  used  for  full  scale  predictions  prior  to 
trials. 


Notional  time  histories  of  full  scale  ship  speed  Us{t,)  and  propeller  RPM  Ns(ts)  as  a 
function  of  full  scale  time  ts  are  sketched  in  Fig.  8.  Prior  to  execution  of  a  stemplane  jam 
recovery,  the  steady  ship  speed  and  RPM  during  the  approach  are  denoted  by  t/J()  and  Ns0 .  Let 
ts0  denote  the  time  when  the  RPM  reversal  is  executed,  and  tsi  the  time  when  the  propeller 
speed  drops  to  zero.  At  time  ts2,  the  propeller  reaches  a  new  prescribed  speed  of  Nsrev  (the 

subscript  rev  refers  to  the  reversed  negative  RPM)  and  is  then  maintained  at  this  constant 
negative  value.  As  a  result,  the  ship  experiences  a  sharp  drop  in  ship  speed  as  noted  in  Fig.  8. 


Figure  8.  Full  Scale  Ship  Speed  and  Propeller  RPM  Reversal 
Given  Us(ts)  and  Ns(ts) 

Existing  RCM  Test  Procedure  for  Jam  Recovery  Simulation 

As  in  a  full  scale  jam  recovery,  a  steady  approach  speed  is  required  for  RCM  correlation 
maneuvers.  The  scaled  model  approach  speed  must  be  set  at 
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(33) 


U 


mO 


"7T 


where  X  =  LjLm  is  the  length  scale  ratio  between  full  and  model  scales,  and  we  will  defer  the 
explanation  of  Eq.  33  until  a  later  section  in  Eq.  43. 

Now,  let  J  denote  a  propeller  advance  coefficient  defined  as  J  =  U/NDp ,  where  Dp 
represents  the  propeller  diameter.  In  an  ideal  fluid,  the  propeller  RPM  required  to  match 
advance  coefficients  at  each  scale, 


u 

'  U 

mp 

m 

NDp 

(34) 


is  then  expected  to  produce  the  same  thrust  coefficient  at  each  scale.  However,  in  a  real  fluid,  as 
a  result  of  viscous  effects,  the  ship  hull  drag  coefficient  Cr  is  greater  at  model  scale  than  it  is  at 

full  scale,  ( [Cr]m  >  [Cf]  J.  So,  the  propeller  must  produce  additional  thrust  to  overcome  the 

additional  drag  at  model  scale.  In  other  words,  the  propeller  requires  additional  RPM  in  a  real 
fluid  than  for  the  ideal  fluid  case,  and  therefore  must  operate  at  a  lower  propeller  advance 
coefficient,  Jm<Js- 

The  required  RPM  at  model  scale  is  given  by 

N„  =  K0Nj I,  (35) 

where  K0  is  a  coefficient  such  that  K0  >  1 .  The  value  of  K0  depends  on  the  hull  form  geometry 
and  ship  speed.  Typically,  K0  is  obtained  from  model  tests. 

At  time  tm0  the  propeller  reversal  is  executed.  (Note  that  the  relationship  between  model 
scale  time  tm  and  full  scale  time  ts  will  be  provided  in  a  later  section.)  At  time  tm2  and 
thereafter,  propeller  rotation  speeds  in  the  existing  procedure  are  set  by  using 

Nmnv  =  NsnvyfI,  (36) 

Shown  in  Fig.  9  is  a  time  history  of  model  speed  resulting  from  the  above  procedure  that 
has  been  (Froude)  scaled  to  the  full  scale  range  and  then  compared  with  the  target  full  scale  time 
history  to  be  matched.  As  noted  in  the  sketch,  the  time  history  of  predicted  ship  speed  can  differ 
from  the  target  ship  speed  by  as  much  as  12%  to  14%.  Since  dynamic  forces  are  proportional  to 
speed  squared,  significant  error  in  the  simulation  of  hydrodynamic  forces  and  moments  can  be 
present  in  RCM  correlation  tests.  This  problem  results  from  the  fact  that,  while  the  approach 
speed  is  properly  Froude  scaled,  the  necessary  Froude-scaled  speed  at  each  instant  of  time 
throughout  the  maneuver  is  not  preserved  during  RCM  tests.  A  method  to  preserve  the  Froude- 
scaled  speed  during  RCM  testing  is  presented  next. 
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12-14% 


Figure  9.  Model  Measured  vs.  Target  Ship  Speed 


Proposed  RCM  Test  Procedure  for  Jam  Recovery  Simulation 

The  procedure  to  establish  the  approach  speed  is  the  same  for  the  new  method  as  it  is  in 
the  existing  approach.  However,  the  reversal  procedure  is  different.  Consider  the  motion  along 
the  body  x-axis.  Let  T  and  t  denote  propeller  thrust  and  thrust  deduction.  A  negative  sign  is 
assigned  to  the  reversed  thrust  to  indicate  that  it  now  acts  in  the  same  direction  as  the  body 
resistance  since  the  vehicle  is  still  moving  forward.  Since,  the  thrust  deduction  is  not  sensitive  to 
viscous  effects,  we  find  that 


-7X1-0-/? 

\/2pU2S 

Expressed  in  dimensionless  form,  we  obtain 

(38) 

During  the  early  stage  of  the  RPM  reversal  while  the  vehicle  still  maintains  an  appreciable 
forward  speed,  [cjm  >  [Cff]j  and  1  -  /  >  0 ,  which  then  dictates  that 


-  7(1-0-/? 

m 

\/2pU2S 

(37) 


[QL>[CrL-  (39) 

So,  we  find  that  the  reverse  thrust  produced  by  the  model  is  too  large,  and  as  a  result  of 
scale  effects  acting  to  increase  model  resistance,  less  negative  thrust  is  required  during  the  early 
stage  of  the  RPM  reversal.  Instead,  we  propose  to  modify  Eq.  36  with 

Hmrev  =  K(ts)Nsrevy[A,  (40) 

where  K(ts)  is  a  time  varying  coefficient. 
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The  RCM  model  is  operated  under  digital  closed-loop  control  with  two-way 
communications  from  the  vehicle  to  an  operations  center  accomplished  by  using  radio  signals. 
Therefore,  the  ship  speed  and  RPM  data  can  be  up  linked  in  real  time.  This  means  that  RPM  can 
be  adjusted  during  a  maneuver  to  ensure  that  the  Froude-scaled  model  speed  agrees  with  the 
target  speed.  The  behavior  of  K(ts )  is  likely  to  depend  on  hull  form  and  Reynolds  numbers,  and 

validation  with  experiments  will  be  required  to  determine  the  shape  of  the  K(ts)  function. 

During  the  later  stage  of  the  RPM  reversal,  ship  speed  has  decreased  from  the  level 
maintained  during  the  approach.  Furthermore,  the  negative  reversed  thrust  produced  by  the 
propeller  will  now  be  substantially  greater  than  the  ship  resistance.  Hence,  the  needed  thrust 
coefficient  is  likely  to  now  be  the  same  at  model  and  full  scale.  We  expect,  therefore,  that 
— >  1  towards  the  end  of  the  maneuver.  A  sketch  of  the  modified  shape  of  the  RPM 

function  is  provided  in  Fig.  10. 


Targ*.  «-L»M  s 


Simulated: 


Figure  10.  Proposed  Procedure  to  Set  Model  RPM 
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Summarizing,  the  existing  approach  for  establishing  the  RPM  level  in  crashback 
simulation  is  based  on  using  the  same  propeller  advance  coefficients  for  model  and  full  scale 
propellers,  namely  Jm  =  Js .  In  the  proposed  approach,  the  RPM  will  be  adjusted  such  that 

Froude-scaled  speed  is  preserved  during  the  maneuver.  We  will  see  in  the  next  section  that  this 
also  has  the  benefit  of  preserving  the  Strouhal  scale  as  well. 

Transfer  Functions  (Scaling  Formula) 

The  fundamental  purpose  of  model  testing  is  to  use  the  model  data  to  predict  full  scale 
performance.  The  formulas  to  relate  model  and  full  scale  ship  performance  are  termed  transfer 
functions  or  scaling  formulas.  These  are  now  derived. 

The  dynamic  equations  for  linear  and  angular  motions  are  given  in  Eqs.  2.  We  postulate 
that  if  the  LHS  of  equations  representing  the  forces  and  moments  acting  on  the  ship  are 
simulated  in  model  tests,  namely  if  LHSm  =  LHSs ,  then  the  motions  experienced  by  the  ship  will 

be  duplicated  by  the  model  RHSm  =  RHSs .  The  relationship  between  full  scale  and  model  scale 

behavior  can  be  deduced  from  these  response  functions.  The  derivation  of  the  transfer  functions 
will  start  with  linear  motion. 


Linear  Acceleration  Scaling 

Consider  the  jc-component  of  linear  motion  given  in  Eqs.  4 


/  \d2 

K  +  roJ^2 

(  »2  \ 
/  \d  x 

K  +  w,)^2 

1/2  pU2S 

1/2  pU2S 

V  ) 

m 

K  j 

(4) 


Let  a  =  d2x/ dt 2  denote  the  linear  acceleration.  Although  added  mass  is  associated  with  the 
transformation  of  kinetic  energy  in  ship  motion,  viscous  effects  associated  with  this  quantity  are 
small  [Szebehely,  1952,  1953].  Equation  4  may  be  simplified  to 


as  _  K  +  ^L 

pU2S 

s  KL 

K4  .  u;  l. 

K+«,L 

pU2S 

«  KL 

\pU2s]m  U2m  Ls 

(41) 


A  comparison  of  Eqs.  41  and  27  suggests  that  the  solution  for  linear  acceleration  can  be  set  by 

-344  =  1  and  U2S=U2±.  (42) 

am\L)  4, 

We  see  that  the  time  history  of  linear  acceleration  measured  at  model  scale  is  directly  applicable 
to  full  scale  with  the  proper  time  scale  applied. 
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Linear  Velocity  Scaling 

To  satisfy  Eqs.  27,  32  and  42  requires  that  the  model  be  tested  at  a  speed  given  by 

=  £/.(/.),  (43) 

where  A  =  LjLm  is  the  length  scale  ratio  between  full  and  model  scales,  and  Um(tm)  is  termed 

the  target  speed.  Note  that  this  relationship  must  be  satisfied  throughout  the  entire  time  history 
of  the  maneuver.  The  Reynolds  scale  effect,  which  alters  model  scale  resistance,  requires  that 
propeller  RPM  be  adjusted  to  compensate  and  preserve  this  Froude-scaled  speed  relationship  as 
discussed  in  the  previous  section.  Note  also  that  Eq.  43  is  identical  to  Eq.  32,  which  is  one  of  the 
two  conditions  to  also  satisfy  the  Strouhal  number.  The  other  condition,  Eq.  31,  will  be 
addressed  below  in  the  section  on  angular  velocity  scaling. 

Time  Scaling 

By  definition, 


dts 


(44) 


Integrating,  we  find  that 

\dts=\-dUs 

J  J°s 

ts  -  [—  VA  dUm  where  as  =  am  and  d Us  -  -Jx  dUm ,  (45) 

iam 

t  —  yfX  t 

s  m 

and  this  represents  the  relationship  between  the  time  base  to  be  used  at  full  and  model  scales. 

Linear  Displacement  Scaling 

In  a  similar  manner,  we  have,  by  definition, 


dt. 


(46) 


Integrating  this  expression,  yields 

\dxs  =  \U,  dts 

Xs(t, )  =  \rx  Umy[I  dtm  where  Us=y[XUm  and  ts=s[Xtm,  (47) 

xs(ts)=  Xxm{tm) 
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and  relates  the  time  history  of  ship  trajectory  at  each  scale. 


Angular  Acceleration  Scaling 

Recall  from  Eqs.  4,. 
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As  with  added  mass,  we  will  assume  that  viscous  effects  associated  with  added  moment  of 
inertia  are  small,  that  is 


[/0  +  [\  ]  m  _  [A)]m 

['.+/,],  “FT 


(48) 


Moment  of  inertia  I  is  defined  by 

Ix  =  j x2dm  Iy  =  jy2dm  Iz  =  jz  dm  . 


We  will  make  the  assumption  that  the  model  can  be  fabricated  to  satisfy 


\j X 1  m  _  ^  ('Jiti  _  [^i  1  m  Ps  ^5 

TErETET"* 


(49) 


(50) 


In  practice,  to  satisfy  Eq.  50  is  not  an  easy  task.  Information  on  full  scale  moment  of 
inertia  data  will  be  required.  Furthermore,  obtaining  a  required  mass  distribution  in  the  model 
may  prove  to  be  time  consuming  or  impossible  to  achieve.  A  method  to  circumvent  this  possible 
fabrication  difficulty  is  being  investigated. 

Substituting  Eqs.  48  and  50  into  Eq.  4  gives 


d2a 

1 

d2a 

dt 2 

"  A 

s 

dt2 

(51) 


This  equation  states  that  the  angular  acceleration  time  history  observed  at  full  scale  can  be 
predicted  from  model  angular  acceleration  data  divided  by  the  scale  ratio  A  =  LjLm  . 

Angular  Velocity  Scaling 

The  angular  velocity  is  defined  as 
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da  ed2a  , 
—  =  — rdt. 
dt  J  dt 2 


(52) 


Note  that  dts  =  -JJ,  dtm  is  invariant  in  linear  motion  or  in  angular  motion.  Substitution  of  this 
and  Eq.  5 1  into  the  above  equation  results  in 


da 

1 

da 

.dt . 

^dt _ 

(53) 


This  is  an  important  result.  One  of  the  assumptions  required  for  the  Strouhal  scale  to  be 
preserved  during  model  testing  is  Eq.  31,  which  is  repeated  below. 


da 

da 

.dt . 

m 

.dt  _ 

(31) 


Substituting  Eq.  43  into  Eq.  3 1  shows  that  it  matches  Eq.  53.  Therefore,  if  these  relationships 
obtained  from  Froude  scaling  are  preserved  throughout  the  maneuver,  then  the  Strouhal  scale 
will  also  be  matched  ensuring  that  unsteady  effects  are  similar. 


Angle  Scaling 

An  angle  results  from  integration  of  an  angular  velocity  as  follows 


Substitution  of  dts  =  -Ja  dtm  and  Eq.  53  into  the  above  equation  gives 


(55) 

This  is  an  important  result.  One  of  the  assumptions  required  for  restoring  moment  to  be 
preserved  in  RCM  tests  is  the  second  of  Eqs.  27.  This  assumption  is  now  verified  by  Eq.  55, 
where  ai  =  {(p,0,y/)  and  (p ,  6  and  ip  represent  roll,  pitch,  and  yaw  angles,  respectively. 


Summary  of  Transfer  Functions 

The  scaling  formulas  that  relate  model  and  full  scale  maneuvering  behavior  are 
summarized  in  Table  3.  We  are  pleased  to  report  that  these  formulas  are  used  in  the  existing  test 
procedure.  We  have  systematically  derived  these  formulas  from  first  principles  using  the 
similarity  scaling  approach  defined  herein.  The  requirements  needed  to  obtain  Table  3  are 
clearly  specified.  A  new  test  procedure  is  proposed  so  that  these  requirements  are  realized  in 
RCM  tests. 
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Table  3.  Scaling  Formulas 


Quantity 

Formula 

Linear  Acceleration 

Linear  Velocity 

UM=Sum(i.) 

Linear  Displacement 

5 

II 

** 

Angular  Acceleration 

i_ i 

1  ’ 

"a 

J  s  L 

1  1 

M  42 

m 

Angular  Velocity 

da 
dt  _ 

1 

,~VI 

da 

~dt_ 

m 

Angle 

Time 

ts=^L 

Within  the  framework  of  the  present  mathematical  formulation,  we  have  shown  that  the 
Strouhal  scale  can  be  satisfied  in  RCM  tests  if  the  Froude  scale  is  matched.  Corrections  have 
been  introduced  in  order  that  Reynolds  scale  effects  on  forces  and  moments  acting  on  the  full 
scale  craft  are  correctly  simulated  in  model  tests  with  the  proviso  that  the  assumptions  introduced 
in  Eqs.  48  and  50  hold.  To  compensate  for  the  scale  effect  on  ship  resistance,  propeller  RPM  can 
be  adjusted  to  match  the  Froude-scaled  target  speed.  Furthermore,  as  a  result  of  relatively  low 
Reynolds  numbers  on  the  model,  and  the  location  of  appendages  in  a  relatively  thick  model  hull 
boundary  layer,  significant  scale  effects  on  forces  and  moments  will  be  experienced  by:  rudders 
and  stemplanes,  the  sail  (moderate  effect),  and  the  ship  hull  (smaller  effect).  Methods  to 
alleviate  these  effects  on  forces  and  moments  have  been  proposed. 

Concluding  Remarks 

We  consider  the  modeling  of  maneuvering  histories  of  ship  motion,  trajectory  and 
attitude  to  be  expressed  by  differential  equations  with  terms  that  resemble  a  spring-mass  system. 
Forces  and  moments  acting  on  the  vehicle  and  the  subsequent  ship  motions  in  RCM  tests  and  full 
scale  trials  are  formulated  based  on  a  similarity  approach.  Within  the  framework  of  the  present 
mathematical  formulation,  requirements  and  test  methods  to  match  Froude,  Strouhal,  and 
Reynolds  scales  (achieved  dynamically)  in  RCM  tests  are  discussed. 

We  postulate  that:  if  the  Strouhal  scale  can  be  preserved,  then  the  unsteady  forces  and 
phase  produced  by  the  ship  in  response  to  the  excitation  forces  and  the  subsequent  unsteady 
motions  will  be  correctly  simulated.  Studies  reveal  that  the  Strouhal  scale  can  be  satisfied  in 
RCM  tests  if:  (1)  the  Froude  scale  can  be  matched  and  (2)  Reynolds  scale  effects  on  forces  and 
moments  acting  on  the  craft  can  be  correctly  simulated  and  (3)  restoring  moments  in  pitch  and 
roll  motions  can  be  correctly  simulated. 
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The  Froude  scale  at  each  instant  of  the  maneuver  can  be  satisfied  by  adjusting  the 
propeller  RPM.  Experience  must  be  gained  to  develop  a  workable  procedure  to  tune  the  RPM 
time  history. 

As  a  result  of  relatively  low  Reynolds  numbers  on  the  model,  and  the  location  of 
appendages  in  a  relatively  thick  model  hull  boundary  layer,  significant  scale  effects  on  forces 
and  moments  will  be  experienced  by:  rudders  and  stemplanes,  the  sail  (moderate  effect),  and  the 
ship  hull  (smaller  effect).  Methods  to  alleviate  these  effects  on  forces  and  moments  have  been 
proposed.  During  a  tight  turn  maneuver,  scale  effects  on  flow  separation  will  have  to  be 
investigated. 

Scaling  formulas  have  been  derived  from  similarity  solutions.  We  have  assumed  that 


Pm 


(50) 


is  satisfied  on  RCM  models.  If  not,  the  scaling  formulas  shown  in  Table  3  require  modification. 
Finally,  the  assumption  that  added  mass  and  moment  of  inertia  are  insensitive  to  Reynolds  scale 
effects  should  be  verified. 
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